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Graphene nanoribbons provide an opportunity to integrate phase- coherent transport phenomena 
with nanoelectromechanical systems (NEMS). Due to the strain induced by a deflection in a graphene 
nanoribbon resonator, coherent electron transport and mechanical deformations couple. As the 
electrons in graphene have a Fermi wavelength \f ~ ao ~ 1.4 A, this coupling can be used for 
sensitive displacement detection in both armchair and zigzag graphene nanoribbon NEMS. Here it 
is shown that for ordered as well as disordered ribbon systems of length L, a strain e ~ (w/L) 2 due 
to a deflection w leads to a relative change in conductance 5G/G ~ (w 2 /aoL). 



Nanoelectromechanical (NEM) resonators hold 
promise for technological implementations such as tun- 
able RF-filters and ultrasensitive mass-sensing. NEMS 
are also of interest in connection with fundamental 
studies of quantum properties of macroscopic systems. 
Regardless of application area, transduction mechanisms 
for system control and readout must be implemented. 

Being only a single atomic layer thick, graphene consti- 
tutes the ultimate material for 2D-NEMS,and graphene 
NEMS have already been demonstrated [l|-|6|. Because 
electron transport through mesoscopic graphene devices 
can be phase coherent [7H9|, using graphene in NEMS 
means that phase coherent transport phenomena can be 
directly integrated into NEM resonators. This allows the 
motion of the NEMS to couple to the length scale set by 
the Fermi wave length ~ &o ~ 1.4 A. 

So far graphene NEMS have operated in the diffusive 
transport regime where electric [5|, |6[ transduction has 
been based on charge carrier density modulation. In 
those experiments the graphene was suspended above a 
backgate a distance d as shown in Fig. QJa). For a sheet 
of length L and width W the capacitance to the gate 
is Cq ~ eoLW/d pH [ll[. Hence, the backgate voltage 
Vq induces a carrier density n^e = epVo/d. This leads 
to a conductivity of a = {leoVo/d [l2|, where fi is the 
mobility. Motion detection then uses the change in car- 
rier density with distance. For a deflection w away from 
the equilibrium distance d, the relative change in con- 
ductance is Serf a ~ w/d [13|- Note that only geometric 
length scales (w and d) enter into this expression. 

A deflection w also induces strain e ~ (w/L) 2 , which 
affects both the dynamical 0, Q and the electronic [Hj 
properties. For diffusive transport, strain leads to a lin- 
ear increase in resistance [16[ with a constant of propor- 
tionality of order unity. Hence, the relative change in 
conductivity due to strain, Sa str ain/o' ~ (w/L) 2 , is typi- 
cally negligible compared to that from the carrier density 
modulation (for a more detailed anaysis, see Ref. [HI). 

The situation is different if one considers coherent 
transport in graphene nanoribbons. The transverse con- 
finement then give rise to conductance quantization. 
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FIG. 1: (Color online) Suspended graphene NEMS systems, 
(a) Graphene sheet of length L and width W suspended in the 
xy— plane above a backgate. (b) Single suspended nanorib- 
bon. For armchair nanoribbons, the strain induces a trans- 
port gap which can be used for displacement detection, (c) 
Graphene nanoribbon interferometer for displacement detec- 
tion using zigzag nanoribbons. 



This prevents conductance changes due to motion in the 
backgate electrostatic field leaving strain as the only cou- 
pling between deformation and conductance. In this pa- 
per it is shown that in such graphene nanoribbon NEM- 
de vices, an operating point where a changes with dis- 
placement as 5 a /a = Sa str ain/o' ~ (w 2 /clqL) can be 
found. For armchair nanoribbons [see fig. DJb)], this is 
due to the opening of the transport gap. For zig-zag 
nanoribbons, which has no transport gap, an interferom- 
eter type set-up [see Fig. DJc)] can be used. This set-up 
can also be used in the presence of edge-disorder. 

Typically, graphene NEMS in equilibrium is not under 
zero strain. Not only will this inhibit ripple formation, it 
will also lead to more linear mechanical response. This 



2 



strain can either be due built in strain or due to biasing 
to a working point wq. In the latter case, the sensitivity 
to a variation w — wq+5w in deflection is naturally linear 
in 5w, i.e. SG/G ~ (w 2 /a Q L) « (wq/L)(5w/clo). 

Transport through suspended graphene sheets and rib- 
bons and in graphene with strained regions has been 
studied previously by several researchers [171-120] and the 
prospects of using strain in a controlled way to influence 
electronic properties is currently an active research field. 
Here the focus is on displacement detection in graphene 
nanoribbon NEMS. 

The electronic properties close to the charge neutrality 
point are well described by the nearest neighbor tight 
binding model. Suppressing spin indices it is 

n,5i rn 

(i) 

Here is the creation operator for an electron at the 
point R n = (niai + n 2 a 2 ) and b n s i the destruction op- 
erator for electrons at the site R n + Si. The basis is 
here: ai = a (3/2, x/3/2), a 2 = ^(3/2,-^3/2), S 1 = 
a (l/2, y/3/2), S 2 = a (l/2, -y/3/2), and S 3 = a (-l, 0). 
For unstrained graphene, t n s = to « 2.7 eV, and the 
Fermi velocity is Tivp = 3£o&o/2. 

For uniform strain e along the x— direction (armchair 
edge) the bond-lengths change from ao to 

N = \S 2 \ = ao[l + 0.256(1 - 3<7 P )], |5 3 | = a [l + e], 

while for uniform strain in the ^-direction (zig-zag edge), 

\Si\ = \S 2 \ = ao[l + 0.25e(3-(7 p )], \S 3 \ = a [l - ea p }. 

Here <j p « 0.1 is the Poisson ratio. 

The changed lengths alter the hopping energies as 
t n s i = to(l + A^). Typically, e <C 1 and it suffices to 
work to first order in e. As A^ oc e, only first order terms 
in Ai need to be kept. The spectrum then remains gap- 
less and linear and can, for uniform strain, be described 
by the low energy Hamiltonian H = Hp + V where 
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Here E is a modified a-matrix defined as 
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k = ^-(V3[A 2 -Ai], Ai+A 2 -2A 3 ). 



Hence, the Fermi velocity changes and becomes 
anisotropic, and the locations of the Fermi points in 
wavevector space changes. 

The A^:s depend on the direction the strain is ap- 
plied in and must be determined from first principles. 
In the context of carbon nanotubes, this has been stud- 
ied extensively 2ll-[23| mainly using tightbinding Hiickel 



theory or Koster-Slater calculations [24|. More recently, 
density functional theory has been applied to strained 
graphene H, [26j. Here, the model of Ribeiro et al. (26| 
will be used, where ts i = to exp[— j3i(Si/ao — 1)]. 

For strain along the x— direction (armchair) /3i = /3 2 = 
2.6, /?3 = 3.3 [2^], and to lowest order in e, one finds 

E = (1 - 2.3be)a x x + (1 - 0A6e)a y y, (3) 
2 

ko = 2.84q/. (armchair) (4) 

3ao 

Here a XiV are the conventional Pauli spin-1/2 matrices. 
For strain along the y— direction (zig-zag) j3\ = f3 2 = 3.15 
and /3s = 4 [26[ , which gives 



E = (1 - 0.74e)a x x + (1 - 2.30e)a y y, 
2 



k 



3ao 



2.70ey. 



(5) 

(zig - zag) (6) 



Consider now an armchair graphene nanoribbon of 
length L, uniform width W, suspended above a back- 
gate as in Fig. QJb). The supported parts (regions I and 
III) are assumed to be unstrained, while the suspended 
part (region II) is under finite strain e > 0. The con- 
ductance in the linear response regime is related to the 
transmission function T{E) as G = 2(e/ti 2 )T(E) 1 where 
the prefactor 2 accounts for spin. 

Confinement in the ^/-direction leads to quantization 
of transverse wavevector components q n — 2i\j3\[3ao + 
nir/W + ko for integer n. Here ko is given by Eq. ((4]). 
If the interfaces between strained and unstrained regions 
are along the ^-direction transverse mode number n will 
be conserved. Hence, T = ^2 n T n and the problem re- 
duces to solving the ID Dirac equation 



(7) 

Here v(x) = V(x)/Hvf is the effective potential in the 
ribbon and £ n = E n /(hvF)- 

To calculate T Eq. ([7]) should be solved in the regions 
I, II and III [see Fig. DJb)] and the solutions matched at 
the interfaces (see also Ref. 1270 . In a region of constant 
v and e, the solution with energy £ in band n is 



ip n (x) = A n e 



ikx 
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B n e 



-ikx 



1 

-i0 n (k) 



.(8) 



a 2 . 



where exp[i<9 n (/c)] = and k = +y/{£ - v) 2 

The matching of wavefunctions between regions is de- 
termined by current conservation [28j. The current- 
operator corresponding to the Hamiltonian in Eq. (0 is 
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FIG. 2: (Color online) (a) Transmission probability T through 
the metallic armchair ribbon (W — 17 nm, L — 28 nm) 
[see Fig. \Hb)] as function of £ and e. The solid lines fol- 
low from the long wavelength approximation [Eq. J9}] for 
the strains e = 0.0% (blue squares) , 0.2% (red circles), 0.5% 
(green triangles), 1.0% (black stars) . The discrete symbols 
were obtained using numerical tightbinding calculations, (b) 
Transmission probability T as function of strain e for ener- 
gies £ = 0.0, 0.01, 0.02, 0.03, 0.04 and 0.05. The slope of 
the dash-dotted line determines the sensitivity of the working 
point around e — 1.0%. 



J x = 2vfYi x o x . Consequently, at the interface between 
regions I and II y // £|V^ = V^V^n However, the fac- 
tors cancel in the final expression for T n 



T n {E) 



sin 



sin V n i - sin U n n 



l 2 N 



COS 6 n i COS 6 n ii 



(9) 



Here (j) n = kuL while # n are the propagation angles 
for electrons in regions I and II. 

In Fig. [2] (a), T is shown for strains e = 0.0 — 1.0% as 
the solid lines. The discrete symbols, were obtained nu- 
merically using the tightbinding Hamiltonian in Eq. (pQ) 
and the relation T = Tr\T lG^T rG^]. Here G^^ are 
the retarded (advanced) Green's functions for the ribbon 
and Tl,r are self energies accounting for semi-infinite 
graphene leads. Also, in the numerical calculation, no 
linearization in strain has been made. As can be seen, 
for the lowest plateau, a transport gap opens up with 
increasing strain. 

From Eq. ([9]), the sensitivity of the conductance G oc 
T to ribbon displacements w can be obtained. For the 
lowest transverse mode q = and for viji = one finds 



To(£) 



£ 2 



k 2 



£ 2 



fcg COS 1 



h 2 



(10) 



This dependence of T on e is shown in Fig.[2fb). Differ- 
ent curves correspond to different back-gate bias points, 
i.e. different values of £ = 0.0,0.01,.... From this fig- 
ure it is clear how for a given strain, one may chose a 
working point £ o (by gating the structure) such that the 
slope of the T(e)-curve is maximal. This maximal slope, 
\dT/de\ max then defines the sensitivity. 
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FIG. 3: (Color online) (a) Transmission T as function of 
energy £ and strain e through the interferometer shown in 
Fig. [He). Broad dark bands correspond to destructive in- 
terference, whereas broad bright bands to unit transmission. 
The figure was generated using the tightbinding Hamiltonian 
and recursive lattice Green's function on a system of width 
W — 2 x 7.5 nm, L — 24 nm, and a 'gap' between the ribbons 
of 0.5 nm. (b) T as function of £ taken along the horizontal 
dashed line in panel a. (c) T as function of e taken along the 
vertical dashed line in panel (a). The thick solid line is the 
result of numerical calculation whereas the thin dashed line 
corresponds to the expression sin0 = sin[1.8e(L/oo)]. 



The smallest sensitivity obtains for the working point 
at £ = [dashed line in Fig.[2jb)]. A lower bound for the 
sensitivity can be found by setting £ = in Eq. ([TO]) and 
solving for the maximum magnitude of the slope. This 
gives |<977<9e| max w ^=(L/a ). Hence, for a deflection 
of magnitude w, the relative change in conductance is 
5G/G = 5T/T~w 2 /(La ). 

This result is valid for a metallic armchair ribbon where 
all edges are perfect and impurities absent. For transport 
restricted to the lowest transverse subband long range im- 
purity scatterers will not affect the transport [29| . How- 
ever, short range potentials will have effect. For armchair 
graphene nanoribbons both theory [30|, and subsequent 
experiments [31( suggest that at low temperature, edge 
disorder induce localization. In this case, transport at 
low energies is goverend by variable range hopping and 
strongly supressed. Hence, schemes relying on a single 
armchair ribbon require nearly perfect edges. 

Zig-zag nanoribbons are less sensitive to edge disor- 
der. However, applying strain will not lead to a trans- 
port gap. Instead, to obtain a sensitivity of w 2 / (Lao) 
an interferometer with the suspended ribbon making up 
one of the arms [see Fig. [He)] can be used. In graphene 
ring-geometries Aharanov-Bohm oscillations have been 
observed at low temperatures [32]. Here, no external B- 
field is required. Instead, the effective gauge field due to 
the strain in the suspended arm is exploited. 

The idea is again to use the lowest quantized conduc- 
tance plateau. For zig-zag nanoribbons this is formed 
from current carried by the edge states. Hence, consider 
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FIG. 4: (Color online) (a) Transmission probability T as func- 
tion of energy £ and strain e through the same structure as in 
Fig. [3l but with added edge disorder (probability to remove 
edge atoms 30%). (b) T as function of £ taken along the 
horizontal dashed line in panel a. (c) T as function of e taken 
along the vertical dashed line in panel (a). The thick solid line 
is the result of numerical calculation whereas the thin dashed 
line corresponds to the expression sin</> = sin[1.8e(L/ao)]. 

an edge-state coming from the unstrained region I, which 
split into the two arms Ha and lib [see Fig. [1] (c)]. The 
state propagating in the strained arm (lib) will aquire 
an extra phase <ft = koL — 1.8e(L/ao) [see Eq. flBJ], and 
consequently one expects interference to modulate the 
conductance with a factor ex sin(0). 

In Fig. [3ja) the result of numerically calculating T 
(using the tight-binding Hamiltonian) from region I to 
region III as function of energy £ and strain e is shown. 
Bright regions correspond to T = 1 and dark regions 
to T = 0. These broad dark regions arise due to de- 
structive interference. The fine structure is the result of 
backscattering at the interfaces where the ribbon split 
(inter- valley scattering). This is also visible in Fig. [3(b) 
where the transmission for a specific strain e = 1.2% is 
shown as function of £ . 

In Fig. [3fc) T is shown for a fixed £ as function of 
strain (thick solid line). The period of the conductance 
oscillations agree well with the plotted function 1 + sin <j) 
(thin black line). Hence, approximating T ~ 1 + sin</> 
yields, as was the case for the armchair ribbon, a sensi- 
tivity to deflections of ST/T ~ (w 2 /Lao). 

The effect of edge disorder on the interference pattern 
is shown in Fig. HJ^a). Here, disorder has been accounted 
for by removing the outermost atoms with probability 
p = 0.3 at random along each of the four zig-zag edges 
of the interferometer. Note that although T as func- 
tion of £ [panel (b)] is highly irregular, T as function 
of strain [Fig. 3] (c)] show clear conductance modula- 
tions. Furthermore, comparing the solid line and the 
dashed line in Fig. [4] (c), it is clear that the the expres- 
sion sin^ = sin[1.8e(L/ao)] is still valid. 

In conclusion, by exploiting the possibility to directly 
integrate coherent electron transport with graphene 



nanoribbon NEMS, the conductance through the struc- 
ture can be made to depend on the mechanical deflection 
w as SG/G ~ (w 2 )/(Lao). This is due to the strain in- 
duced shift of the Fermi-points (synthetic gauge field). It 
is this shift which causes the length scale ao ~ \ F ~ \/kp 
to enter the expression for SG. 
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